Electric field induced nucleation: An alternative pathway to metallic hydrogen 
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Electric field induced nucleation is introduced as a possible mechanism to realize a metallic phase 
of hydrogen. Analytical expressions are derived for the nucleation probabilities of both thermal 
and quantum nucleation in terms of material parameters, temperature, and the applied field. Our 
results show that the insulator-metal transition can be driven by an electric field within a reasonable 
temperature range and at much lower pressures than the current paradigm of P > 400 GPa. Both 
static and oscillating fields are considered and practical implementations are discussed. 
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Metallic hydrogen (MH) was predicted by Wigner and 
Huntington [l| in 1935 and it continues to be on the fore- 
front of scientific inquiry; new predictions include high-T 
superconductivity [2], |3j and novel ordered quantum fluid 
behavior [4j, |5j . A significant body of experimental work 
and advanced numerical simulations [Io| have elu- 
cidated the properties of hydrogen at high pressure and 
have provided a great deal of insight to its rich phase 
diagram llj. From a practical standpoint, if metastable 
MH 12| were to be obtained, it could be a revolutionary 
clean energy source fl3j |. 

Attempts to produce MH have focused on high pres- 
sure techniques. Solid MH has not yet been observed un- 
der static pressures of up to 342 GPa 14, lj|. Dynamic 
compression beyond 200 GPa has also been employed 
d, 0. The only direct evidence so far was the brief 
observation of a highly conductive liquid phase under a 
shockw aye p ressure of 140 GPa and temperatures around 
3000 K 1 1 711 . Recent numerical studies have calculated a 
solid |l8l |20| metal-insulator transition pressure > 400 
GPa and a first order liquid-liquid [19] metal-insulator 
transition pressure ~ 200 GPa near 1000 K. 

In this work, we introduce electric field induced nucle- 
ation (FIN) as a possible alternative means by which to 
realize MH; both thermal and quantum FIN are consid- 
ered. The underlying thermodynamic theory is general 
in nature and not concerned with the microscopic aspects 
of the phase transition, whether it occurs by dissociation 
of the molecular form to a monatomic metal, as origi- 
nally considered [l[ , or from a proton- paired insulator to 
a molecular metallic state [2(j. Beyond MH, FIN could 
open new venues for phase transformations to otherwise 
unobtainable materials. 

FIN is a recently developed concept of metal phase 
nucleation in an insulating host under a strong static 
21 1 or oscillating [22| field. FIN produces needle-shaped 
conductive embryos with vertical axes aligned with the 
field (see Fig. [IJ. It has been identified as a mechanism 



for crystallization in chalcogenide glasses [2l|. Also, it is 
supported by observations of laser induced nucleation of 
particles aligned with the beam polarization 23] . 

Given the high electrical conductivity of a metallic nu- 




FIG. 1: Electric potential surface plot and field streamlines 
show the greater electrostatic energy reducing effect of elon- 
gated metallic nuclei versus spheres. The corresponding de- 
crease in the the nucleation barrier W/Wo with respect to the 
applied field E/Eq is shown in the lower plot. 



cleus, the field E causes it to act as a strong dipole, 
p cx E. The associated negative contribution, ~pE, to 
the free energy suppresses the nucleation barrier. The 
premise of FIN is that anisotropic particles can result in 
much greater barrier suppression than spheres, as shown 
in Fig. [1] Specifically, a metal particle of volume f2 in a 
uniform field E reduces the free energy according to [2J] , 
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where e is the permittivity of the insulating host and n 
is the depolarizing factor, which depends on the particle 
geometry (see Fig. [I}. For example, a sphere provides 
n = 1/3, while a prolate spheroid of radius R and height 
H yields n = (R/ H) 2 [ln(2H/ R) - 1]. An approximate 
relation n ~ (R/H) 2 holds for arbitrary needle-shaped 
particles. Indeed, the condition that the dipole field bal- 
ances the external field in a metal determines p ~ EH 3 , 



2 



leading to F E pE ~ E 2 H 3 &E 2 {H/R) 2 with 

n ~ r 2 h. 

According to FIN, the free energy of a MH nucleus can 
be expressed as, 



F = A<t + Vl\h\+F e , 



(2) 



where A is its surface area, a is the surface tension, and 
p is the chemical potential difference between the two 
phases. The positive chemical potential term in Eq. ([2]) 
corresponds to a system that is stable in zero field. 

The exact shape of the elongated nuclei is not known, 
but modeling with either spheroidal or cylindrical par- 
ticles leads to differences only in numerical coefficients 
21[. We opt for the mathematically more concise form 
of a cylindrical nucleus with the free energy of Eq. ([2]) 
given by, 
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Here, we have approximated the logarithmic term in Eq. 
(jTj) with unity, and introduced the nucleation barrier, 
Wo = 167T(7 3 /(3p 2 ), and radius, Ro = 2a/\fj,\, for spheri- 
cal particles of the classical nucleation theory; the char- 
acteristic field is E Q = 2^W /sR? ) . In Fig. El a contour 
plot of the free energy illustrates how the system can 
lower its free energy by forming elongated particles. 
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FIG. 2: (a) Contours of free energy from Eq. positive and 
negative regions separated by the zero contour. Nucleation of 
elongated particles of critical height H c /Ro is governed by the 
barrier, W. (b) Cross-section of the free energy from (a) along 
the line a shown with (solid) and without (dashed) the electric 
field E/Eq. Thermal and quantum nucleation pathways are 
labeled Pt and Pq. 

The free energy of Eq. ^ seems to suggest that nuclei 
with R — > are energetically favorable. Realistically, R 
must be greater than some minimum value determined by 
extraneous requirements, such as sufficient conductivity 
to support a large dipole energy or mechanical integrity. 
Based on data for other types of systems, it was estimated 
2l| that a reasonable minimum radius is R m in = aRo, 
where a ~ 0.1. Lacking other data we maintain a ~ 0.1, 
which leads to R m in on the order of interatomic distance. 
The free energy in the region R/ Ro < a is substantially 



larger than described by Eq. (|3|), since the energy re- 
ducing effect of the electric field cannot be manifested 
by such small particles. As a simplifying approximation 
we consider nucleation along the path R/Rq = a (see 
Fig. E]); alternative paths that start from the point (0,0) 
introduce only insignificant numerical factors. 
From Eq. |3]), the nucleation barrier is, 



W = (E C /E)W , E c 



* 3/2 E . 



(4) 



Thermally activated FIN, with the probability Pt oc 
exp(—W/kT), becomes important when E > E c = 
c?l 2 E§. The requirement of an elongated nucleus, 
H/R ^> 1, places an upper limit on the field, E < E c /a 2 . 

For MH, using ^ [i ~ 0.1 Ry/atom - 10 11 J/m 3 and 
a ~ 1 Ry/atom (for structurally different phases) yields 
Ro ~ 10 A and W - /i#o ~ !0 3 e V. Then with a = 0.1, 
and [9(] e ~ 1, we obtain E c ~ 10 7 V/cm. 

Following the approaches in Refs. 0, 27, quantum 
nucleation has the probability exponent, 



2M(R,H)F cyl (R,H)ds, 



(5) 



where M(R,H) is the local mass. The integral is taken 
over the tunneling path C with element ds through the 
two-dimensional free energy barrier. 

To determine M, we consider a host material of density 
p in which a cylindrical nucleus of density p + Ap grows 
by the inward flux of particles with radial speed, 

v r = - — -R for H^>r>R, (6) 
p r 

where R is the radial growth velocity of the nucleus. The 
associated hydrodynamic kinetic energy takes the form, 



H 



E k a (p/2) / v 2 (2irH) rdr = MR 2 /2, 
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from which, 



M(R, H) = 2nR 2 H— In ( — 
p \R 



(7) 



(8) 



Here, we have taken into account that the assumption 
of an axially symmetric velocity field in Eq. (JBJ) is valid 
when r <C H . The region r 3> H is closer to spherical 
symmetry (considered in Refs. [26|, [27[) where the inte- 
gral for kinetic energy is dominated by its lower limit and 
is relatively insignificant. 

Inserting Eqs. and ([8} into Eq. (J5]) and using 
dimensionless variables x = R/Rq and y = H/Rq, results 
in the barrier factor, 
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ds, (9) 
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where So = (6irAp 2 R^W / ph 2 ) 1 / 2 and the integral is 
dimensionless. Optimizing the latter leads to an equa- 
tion of motion for the optimal path that is analytically 
intractable. We limit our approximation to the suffi- 
ciently simple but physically reasonable path along the 
line x = a -C 1, as shown in Fig. [5] The path ends 
at the point Ho/Ro = V3E c /(aE), which follows from 
F cy i{Rmin,Ho) = 0. Evaluation of Eq. ((9|) then yields 
the quantum nucleation probability Pq ex. exp (-S), 
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where 7 = [a/ In (V3E c /a 2 E)] ' *. The field-dependent 
expression in parentheses describes the FIN effect; setting 
that expression to unity reduces Eq. (fTU]) to that derived 
in Ref. [26| (to within insignificant numerical factors). 
Several other pathways through the barrier differing from 
the line x — a were found to result in larger values for the 
barrier factor, making our approximation along the line 
x = a a reasonable estimate of the variational minimum. 

Neglecting any difference in the probability pre- 
exponentials, the temperature of interplay between ac- 
tivated and quantum nucleation is given by 



Wo l 2 E 
kSo E c 



fll) 



If, for example, we consider the solid molecular to solid 
metallic transition near T = 0, then [l| Ap/p ~ 10 and 
Ap ~ 1 g/cm 3 . With the previous values of Wo ~ 10 3 
eV and Ro ~ 10 A, the temperature regime of quantum 
nucleation is shown in Fig. Eta). 
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FIG. 3: (a) Temperature regime of quantum nucleation (T < 
Tq) as a function of applied field E / E c . (b) An approximate 
T — E phase diagram of field induced nucleation (with the 
field given by the relative laser intesity E/E c = y/I/I c ) . 

Our treatment thus far equally applies to the case of os- 
cillating laser fields with the equivalent static field given 
by E = \J 4mT/c, where / is the beam intensity and c 
the speed of light. As detailed in Ref. [22J, the equivalent 
static field can be applied under the following conditions: 
1) the radius of the metal particle less than the field de- 
cay length, R < S = c^l/2irauj, where ui is the laser 
frequency and a is the MH conductivity; and 2) the MH 



particle must be sufficiently polarizable to ensure a strong 
field effect. The latter condition can be expressed [22| . 
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where uj p is the MH plasma frequency. With adherence 
to these conditions, all of the above results hold with the 
substitution E/E c — > y 'I /I c where, 
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(13) 



Optical absorption can facilitate a temperature in- 
crease 8T above ambient To so that T = To + ST. In 
terms of material parameters, the activated and quan- 
tum nucleation rates can be estimated respectively as, 



Pt = NqUo exp 



Wn 



k (To + 5T) 



and 



P Q = Af u;oexp ( -S 



Jc_ 

7 2J 



(14) 



(15) 



In the pre-exponential factor No^o 



10 



cm 



No ~ Rq 3 is the number density of fluctuating centers 
and wo ~ Wo/(Soh) is the characteristic vibrational fre- 
quency of such a center. 

For T > Tq, by setting P T = 1 in Eq. (Q3J we ob- 
tain a reasonable transition temperature as a function of 
the applied field (or intensity). The result is an approxi- 
mate T — E phase diagram, as shown in Fig. [31(b), which 
indicates that thermal FIN can be practical. 

Given the very rough numerical values above, we find 
that the probability of quantum nucleation is negligible, 
even at the maximum allowable intensity of / ~ 10 4 / c 
(or static field E = E c /a 2 ~ 10 2 T C ). However, our 
analysis is based on quantum nucleation at T = and 
ambient pressure. At higher temperatures, nucleation 
may be more efficient through an optimum combination 
of activation and tunneling, especially considering the 
broad temperature regime of quantum nucleation [see 
Fig. [3la)]. In addition, since our numerical estimates are 
based on solid molecular to solid metallic transitions at 
atmospheric pressure, phase transitions at higher pres- 
sures and temperatures, such as solid-liquid or liquid- 
liquid, would occur with lower barriers Wo and density 
ratios Ap/p, leading to shorter nucleation times. 

One conceivable implementation of FIN is related to 
recent diamond anvil cell experiments wherein cryogeni- 
cally loaded hydrogen was heated with either continuous 
wave [7[ or pulsed [8( lasers. In those experiments, which 
were designed to investigate the melting line of hydrogen 
up to nearly 150 GPa, a laser of wavelength A ~ 1 /im 
and beam intensity / ~ 10 6 — 10 9 W/cm 2 was focused to 



4 



a spot size of 10-30 firm on a Pt absorber located within 
the hydrogen sample. These melting experiments could 
be slightly modified to assess the FIN theory by increas- 
ing the laser intensity (i.e. increase power and/or reduce 
area) to the range of I > I c while maintaining the fre- 
quency required by Eq. (|T2"1) ; given w p ~4x 10 16 s _1 for 
hydrogen. Metal nucleation could be detected by changes 
in, for example, conductivity, opacity and/or scattering 
within the portion of the sample subtended by the beam. 
Since the electric field is the primary phase change driver, 
the experiments could be conducted at much lower pres- 
sures and over a broad range of temperatures. Simul- 
taneously, one could investigate the correlation between 
nucleus orientation and beam polarization. Similar ex- 
periments could be conducted with hydrogen rich alloys 
such as CH4 (or other paraffins) and SiH4(H 2 )2 [H-31|, 
or LiH„ (n > 1) [13. 

In conclusion, our results show that FIN could present 
a new pathway to metallic hydrogen that is within reach 
of existing experimental techniques. Further investiga- 



tion is imperative since there is theoretical evidence [12 1 
that the metallic phase would exhibit long-term metasta- 
bility upon removal of the field at ambient pressure, 
thereby providing a novel source of clean energy. In addi- 
tion, FIN may be useful in realizing other materials that 
have been thus far unobtainable. 
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